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ABSTRACT 

CXO imaging has shown that equatorial tori, often with polar jets, are very common in 
young pulsar wind nebulae (PWNe). These structures are interesting both for what they reveal 
about the relativistic wind itself and for the (nearly) model-independent measurement of the 
neutron star spin orientation they provide. The later is a particularly valuable probe of pulsar 
emission models and of neutron star physics. We describe here a procedure for fitting simple 3-D 
torus models to the X-ray data which provides robust estimates of the geometric parameters. An 
application to 6 PWN tori gives orientations, PWN shock scales and post-shock wind speeds along 
with statistical errors. We illustrate the use of these data by commenting on the implications for 
kick physics and for high energy beaming models. 

Subject headings: stars: neutron - pulsars - stars: rotation 



1. Introduction 

The Crab nebula has long been known to have 
a subluminous zone and termination shock sur- 
rounding its central pulsar. One of the strik- 
ing successes of the Chandra X-ray Observatory 
(CXO) mission has been to show that this shock 
is likely an equatorial band and that similar struc- 
tures are seen around a number of young pulsars 
(e.g. Weisskopf et al. 2000, Pavlov et al. 2001, 
Helfand, Gotthelf & Halpern 2001, Gotthelf 2001). 
Romani & Ng (2003) argued that the apparent 
symmetry of such PWNe, if interpreted as equa- 
torial tori, allowed a robust fit for the 3-D orien- 
tation of the pulsar spin axis VL and showed that 
measurement of this axis can be effected even for 
quite faint PWNe. They also argued that, taking 
PSR J0538 + 2817 as an example, comparison of 
the spin axis with the proper motion axis v can 
be a sensitive probe of pulsar physics (see Spruit 
& Phinney 1998, Lai et al. 2001). In particular, 
when the alignment is due to rotational averaging 
of the transverse momentum, tight alignments im- 
ply momentum kicks lasting many times the pul- 
sar initial spin period. Such time constraints on 
the kick timescale can exclude otherwise plausible 
models. We wish here to systematize this compar- 



ison by outlining how robust fits for model param- 
eters can be obtained even for relatively low-count 
CXO PWNe images. 

The characteristic scale of the wind termination 
shock around a pulsar of spindown power E is 

r T w (E/incPext) 1 / 2 . 

This structure should be azimuthally symmetric 
about the pulsar spin axis when P ext > P ra m = 
6 x 10~ 10 nw|g/cm/s 2 , i.e. when the pulsar is sub- 
sonic at speed lOOi^km/s. Pulsars will seldom be 
sufficiently slow to show toroidal shocks in the gen- 
eral ISM (where they will have PWN bow shocks), 
but young PSR often satisfy the azimuthally sym- 
metric torus condition in high pressure SNR inte- 
riors, van der Swaluw et al. (2003) have provided 
analytic descriptions of conditions in SNR interi- 
ors, modeling the pulsar parameters required for 
a subsonic PWN as the SNR evolves. 

2. Fitting Model and Technique 

The discovery papers showing obvious toroidal 
PWN structure have, of course, often provided es- 
timates of the torus scale, position angle and in- 
clination. However, most of these estimates were 



simple eyeball matches or, at best 2-D fits of el- 
lipses in the image plane. The ubiquity of clearly 
equatorial structures with toroidal symmetry and 
the robust physical interpretation in terms of a 
static termination shock in an equatorial wind ar- 
gue that direct fits of the implied 3-D structure 
are most appropriate. Further, quantitative com- 
parison with other pulsar observables requires er- 
ror estimates for the geometric parameters, which 
are often not provided. Finally, we seek to con- 
strain the orientations of central tori in PWNe, 
even when they are not immediately obvious. All 
of these considerations require a robust fitting pro- 
cedure with (statistical) error estimates. 

We adopt a default model of a simple equatorial 
torus with a Gaussian emissivity profile in cross- 
section. Clearly, this is not a sufficiently detailed 
model to describe the complete PWN structure in 
the highest resolution images. In fact, for Crab, 
Vela and a few fainter PWNe we find that multi- 
ple tori component significantly improve the model 
fits. Also, we do not fit the detailed radial varia- 
tion of the torus emissivity and wind speed. Such 
extensions of the model could be useful for com- 
parison with numerical models of relativistic shock 
flow (e.g. Komissarov & Lyubarsky 2003, Shibata 
et al. 2003). However, these models are far from 
complete, and for the fainter PWNe the existing 
count statistics do not warrant such model exten- 
sions. Accordingly we adopt the generic torus fit 
and apply this to six pulsars, with increasingly 
poor image statistics. Since for the bright PWNe 
fine structure beyond a simple torus is visible, we 
must assume that these models are incomplete de- 
scriptions of the fainter objects as well. However, 
the simple torus gives an adequate description of 
the data and, assuming that this model captures 
the gross nebular structure, we obtain quantitative 
estimates of the orientation and wind flow param- 
eters. The estimates are quite robust to perturba- 
tions in the input fit assumptions. Note also that 
simpler structures are included as subsets of our 
fit parameter space, e.g. a simple uniform halo, 
which is reproduced by a face-on torus with small 
radius and large blur. Such models are strongly 
statistically excluded, even for the faintest sources. 

We characterize the termination shock (torus) 
by a polar axis at a position angle \I> (0 — 180°, 
measured N through E), with inclination angle £ 
(0 — 180°) from the observer line of sight. The 



torus is simply described by a radius r and a fi- 
nite thickness or 'blur' of Gaussian cross section, 
dispersion 5, about the central ring. We assume 
that the surface brightness (averaged over its cross 
section) is uniform in azimuth. However, the post- 
shock flow is still expected to be mildly relativistic 
and so an azimuthally symmetric ring will vary in 
apparent intensity, due to Doppler boosting. Syn- 
chrotron emission of intensity Iq and photon spec- 
tral index T in the rest frame will have an apparent 
intensity 

Joe (l-n-/3)-( 1+r )/ 

(Polling, et al. 1987) where (3 = v/c is the bulk 
velocity of the post shock flow and n is the unit 
vector to the observer line of sight. 

In practice we set up a coordinate grid with x, 
y oriented with the CCD frame and z along the 
observer line-of-sight. This is rotated to a grid 
aligned with the torus 

x' = —x cos — y sin ^> 

y' = (xsin^ — ycos^) cos£ + zsin£ 

z' = — (xsin ^ — ycos VP) sin£ + z cosC 

with x' along the line of nodes, and z' along the 
torus axis. In this frame the emissivity toward the 
Earth is 

I(x',y',z')= ' 
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where N is normalized to the total number of 
counts in the torus component. We actually as- 
sume a fixed T = 1.5, typical of PWNe and leave 
(3 as a model parameter. 

These parameters define the basic 3-D torus of 
PWN emission. To form the 2-D torus image for 
data fitting, we integrate the emissivity though the 
line of sight to get the counts in each pixel, i.e. 



C{x,y) = Y J I{x',y',z') 



The model image also generally includes a uni- 
form background and a point source for the pulsar, 
whose intensity distribution is the PSF appropri- 
ate to the detector position and source spectrum. 
For Crab and Vela, the point sources have strong 
pile-up and we can optionally blank a small region 
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around the point source. There are a few spe- 
cial cases of this basic model. In several PWNe 
(most notably Vela), the termination shock has a 
double torus, presumably spaced above and be- 
low the spin equator. We model this with identi- 
cal tori, symmetrically offset along the torus axis 
about the pulsar by distance d, following Hclfand, 
Gotthelf & Halpern (2001). In the case of the 
Crab, the structure includes two, apparently inner 
and outer, tori. Finally in several cases relatively 
bright polar jets complicate the analysis. We can 
add model jet emission along the torus axis, blank 
the jet regions or include jet region photons as a 
fixed background. We do not at present fit these 
components separately but their inclusion makes 
for more appealing model images. The torus pa- 
rameters are robust to the inclusion or exclusion 
of the 'jet' region photons. 




Fig. 1. — Torus Geometry - The figure is the plane 
of the sky projection of a double torus with fit pa- 
rameters and angles labeled. £ (not shown) is the 
angle between the sky plane normal (=the Earth 
line-of-sight z) and f2. 

Thus the basic torus model includes parame- 
ters ^, C, r, S, (3 and a background. Additional 
parameters can include the point source counts or 
double torus separation. To keep the model par- 
simonious, we apply the constraint that the to- 
tal model counts Yj xy C xy match the image counts 
(except for the Crab, where bright complex back- 
ground structure is not adequately modeled by a 
constant background). Our fitting procedure uses 
maximum likelihood, with a merit function formed 
from the summed log probabilities of the observed 



counts d xy in each image pixel. We use Poisson 
statistics 



passing to Gaussian statistics for an expected pixel 
count C xy > 20 in the CCD image. In this way 
the distribution of fluctuations in the fitted model 
parameters passes smoothly to \ 2 statistics in the 
limit of high count images. The merit function is 
minimized, using the simulated annealing package 
in Numerical Recipes (Press, et al. 1992). The ini- 
tial parameter estimates are chosen by eye and the 
simplex generally relaxes quickly to the stable final 
solution, with a robust insensitivity to the initial 
conditions. The one exception is the blur thickness 
S, which has a tendency to grow to absorb the un- 
modeled larger scale PWN enhancement present in 
many of the brighter nebulae. As this parameter 
is not particularly useful in the physical descrip- 
tion of the shock, we fixed it at the initial (visual 
estimate) value in a number of cases. The other fit 
parameters are quite insensitive to 5 over a broad 
range 6 < r. 

Equally important are the parameter error es- 
timates. We obtain these via fits to Monte Carlo 
Poisson realizations of the best fit-model. After 
500 random simulations and fits we examine the 
distribution of the fit parameters. These distri- 
butions are well-behaved, with near-Normal dis- 
tributions about the original fit parameters. We 
quote here 1-ct confidence intervals corresponding 
to 68% CL. These are projected (i.e. true multi- 
parameter) estimates of the statistical errors in 
the fit parameters. At very high confidence lev- 
els (> 99%) the parameters become poorly con- 
strained, especially for the fainter PWNe. Never- 
theless we find that the estimated statistical errors 
are quite reliable well beyond 2-er. Systematic er- 
rors, due to jets, complex backgrounds, or other 
unmodeled structures are of course not included, 
but these statistical estimates appear dominant for 
all but the brightest PWNe and thus provide use- 
ful estimates of the parameters of the torus com- 
ponent which captures the gross structure of the 
central PWN. 

The simulations also reveal correlations be- 
tween the fitting parameters. W is, as expected, 
well decoupled from the other parameters. £ and f3 
are substantially co-variant with correlation coef- 
ficients as large as 0.7 since both affect the bright- 
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ness ratio between the front and back sides of the 
tori. Fluctuations toward large inclination thus 
have reduced (3. Similarly for the PWNe with 
large inclinations, £ and r are covariant (coeffi- 
cient 0.3 to 0.7) since the projected angle of the 
torus' bright edge rcos£ is well constrained by the 
images. Our reported la errors include these pro- 
jected correlations. 

Fit parameters and errors are reported for six 
PWNe with a range of brightness in Table 1. We 
start with the 'obvious' tori of Crab and Vela, con- 
tinue with the 'apparent' tori of PSR J2229+6114 
and PSR J1930+1852 and extend to two faint 
'plausible' tori from pulsars for which existing es- 
timates of the proper motion make spin axis mea- 
surements of particular interest. There are per- 
haps 10 additional PWN known where such anal- 
ysis could be fruitful, and we expect there that 
several more useful torus measurements will be 
made as the CXO mission matures. In the present 
analysis, all images are from archival ACIS-S data 
sets, energy cut to 0.5-5kcV. For the Crab we use 
ObsID 1998 (25ks), for Vela we combine ObsIDs 
1987 and 2813-2820 (180ks). For the other objects 
we use the following observations: G54. 1+0.3 (Ob- 
sID 1985, 40ks), PSR J2229+6114 (ObsID 2787, 
lOOks), PSR B1706-44 (ObsID 757, 15ks), PSR 
J0538+2817 (ObsID 2796, 20ks). We mention 
here special features of individual fits. 

The two Crab tori are fitted sequentially, with 
the outer torus fitted first and then used as a fixed 
background for the inner torus fits. The model 
shown combines both components. The fitted tori 
are close to co-planar and have very small sta- 
tistical parameter errors due to the high image 
counts. The Vela tori are spaced symmetrically 
about the pulsar, centered along the symmetry 
axis; we quote the value d in arcseconds projected 
on the sky. For the Crab and for Vela, the high 
S/N images show complicated fine detail in the 
tori and unmodeled fainter structures in the sur- 
rounding PWN. We have tested the fit sensitivity 
to these structures by, eg. blanking the jet regions, 
adjusting the blur thickness and shifting starting 
fit parameters. We find surprisingly small sensi- 
tivity to these changes in the torus fits: the fit 
values are quite robust, although clearly the very 
small statistical errors will be dwarfed even by the 
modest systematic biases. 

For SNR G54. 1+0.3 and PSR B1706-44, the 



global minima of the merit functions pull the tori 
to small radii to absorb excess counts in the PSF 
wings of the central pulsar. By excluding torus 
count contributions to the merit function for pixels 
at < 2", we obtain the minima listed in the table, 
which match the visual structure of the nebulae. 
PSR B1706-44 also has bright regions above and 
below the fit torus, which can be interpreted as 
polar jets and which produce a modest fraction of 
the counts attributable to the putative torus. To 
minimize the sensitivity of the torus fit to these 
'jet' photons, we add these counts as a fixed back- 
ground in the model, as can be seen in the corre- 
sponding model image. 

The PWN of PSR J0538+2817 is quite faint, 
providing only ~ 2% of the point source counts. 
Accordingly, the geometry of the torus is sensi- 
tive to the amplitude of the PSF, and so the point 
source normalization cannot be fitted simultane- 
ously. Instead the PSF amplitude is first fitted 
in the central 1". This (pile-up corrected) PSF is 
held fixed during the torus fitting. Consistency is 
checked by re-fitting for the point source in the 
presence of the torus background, confirming that 
the amplitude is unaffected. Romani & Ng (2003) 
have previously described this PWN, arguing that 
the spatial and spectral distinction from the PSF 
wings make the nebula highly significant, despite 
its low count rate. The fit errors in Table 1 are 
slightly larger that those in Romani & Ng (2003), 
because of the additional fit parameters and the 
multi-dimensional errors. In particular, slightly 
smaller errors (<~ 6° vs. <~ 8°) can be obtained by 
direct fits to ^ alone. However, an equatorial ex- 
tension of the PWN is strongly required in the fit, 
the ^ estimate remains robust to perturbations of 
the starting conditions and a simple uniform halo 
(small £) is strongly excluded - the figure of merit 
degrades to the 3.6a level, referenced to our multi- 
parameter error flags, for models with («0. 

3. Application of Fit Parameters and Con- 
clusions 

In our fitting, we find that (for £ not too small) 
the position angle ^ is the most robust parameter. 
We believe that even if the simple torus model is 
inadequate to fully describe all nebulae, this mea- 
surement of the position angle of transverse PWN 
extension is very robust. Thus comparison of the 
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projected fit axis with the proper motion axis re- 
mains an important application. For Crab, the 
Caraveo & Mignani (1999) HST measurement of 
the Crab proper motion lies at "J/pm — 292 ± 10°. 
This is f 2° or 1.2a off of our inner ring axis. Sim- 
ilarly for Vela, we compare with the new radio 
intcrferometric proper motion of Dodson, et al. 
(2003) at * PM = 302 ±4°, which has both higher 
astrometric accuracy and a corrected treatment of 
Galactic rotation effects from earlier optical esti- 
mates. This vector lies 8.6° from our fitted torus 
axis, a ~ 2.1a discrepancy. Several other compar- 
isons are semi-quantitative at present; we discuss 
these below. 

Neither PSR J1930+1852 in G54. 1+0.3 nor 
PSR J2229+6114 in G106.3+27 has a directly 
measured proper motion. Further, these SNR are 
complex with no clear shell structure, so a velocity 
vector has not been estimated from an offset birth 
site. In contrast, Dodson & Golap (2002) and 
Bock & Gvaramadze (2002) have re-examined the 
controversial association between PSR B1706-44 
and G343. 1-2.3, and argued that this partly shell- 
like SNR is larger than previously believed, placing 
the pulsar well in the interior and making the as- 
sociation more plausible. The former authors find 
a southern extension suggesting a more complete 
circular shell; the vector from the center of this 
structure is at "J/pm ~ 150°. The latter authors 
suggest a proper motion parallel to the expan- 
sion of the nearby bright shell rim, vE'pm ~ 170°. 
These position angles are ~ 10 — 20° from our fit 
axis, but since the SNR evidently suffered asym- 
metric expansion, both of these geometrical esti- 
mates are uncertain. Thus while the axes arc in 
general agreement, a direct measurement of the 
proper motion is essential for any serious compari- 
son. PSR J0538+2817, in contrast, resides in S147 
which has a quite symmetrical structure. Romani 
& Ng (2003) estimated a proper motion axis from 
the offset at -32 ± 4°, i.e. < la off of the PWN 
axis measured here. Kramer, et al. (2003) have 
recently managed to extract a timing proper mo- 
tion for this pulsar which confirms the association 
with S147, although the PA is poorly determined. 
Here, both higher statistics X-ray imaging and an 
astrometric proper motion are needed to effect a 
precision test. 

There are two other young pulsars in shell SNR 
with recent proper motion measurements, PSR 



B1951+32 in CTB80 and PSR B0656+14 in the 
Monogem Ring. Pulsar B1951+32 has a proper 
motion at *ff pm — 252 ± 7° away from its SNR 
birthsitc (Migliazzo et al. 2002) and is presently in- 
teracting with the dense swept-up shell. Thus the 
prominent bow shock seen in HST imaging is not 
unexpected; one would not expect a toroidal wind 
shock. It is plausible that a 'jet' wind could punch 
through the bow shock and Hester (2000) has 
proposed that the Ha 'lobes' bracketing the bow 
shock mark the pulsar polar jets. Under this inter- 
pretation we measure the spin axis is at w 265±5°, 
which is 13° (<~ lAa) away from the proper mo- 
tion axis. 

With a new parallax distance measurement 
Brisken, et al. (2003) and Thorsett, et al. (2003) 
find that PSR B0656+14 is enclosed within the 
~ 66pc-radius Monogem ring. The surprisingly 
small proper motion at ^ pm — 93.1 ±0.4° implies 
a transverse velocity of only 60km/s. If the pulsar 
has a more typical ~ 500km/s space velocity, it 
must be directed along the Earth line-of-sight; in- 
deed, the parallax distance is consistent with the 
near side of the Monogem ring. At its charac- 
teristic age, the pulsar should still be within the 
remnant interior for radial velocities as large as 
~ 600km/s. This SNR exploded in the low density 
local ISM, so we expect the PWN to be toroidal 
with a characteristic radius of ~ 3". Interestingly, 
a short CXO observation shows a faint, nearly cir- 
cular halo around the pulsar at this scale (Pavlov, 
et al. 2002), suggesting a face-on torus. Sched- 
uled CXO observations have the sensitivity to map 
this structure, which we predict will be consistent 
with a torus tilt of <~ 15°. 

Table 2 collects the projected proper motion 
and spin-axis position angle (v, Q) estimates along 
with estimates of the line-of-sight inclination for 
these young pulsars. The trend toward alignment 
(small |A^n.„|) is strong, albeit imperfect. For- 
mally, one should impose a prior on the maximum 
total space velocity v before evaluating the likeli- 
hood of a position angle range {Avp}. If the 3-D 
angle between v and the projected position angle 
is 9, then a physical upper limit on the plausible 
v = uj_cos#/cos\I' restricts the allowed range of 9. 
However, in practice all of these pulsars have rel- 
atively small v — v±, so the disallowed range is 
negligible and the probability of a position angle 
range is simply 2(A^ max - A$ min ) / V . 
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If only the la maximum lA^I are considered, 
then for Crab and Vela alone, there is a 3% chance 
of obtaining alignments as close as those seen from 
isotropically distributed v, Q. However if the 
other three angle estimates of Table 2 are included, 
the chance probability falls to 0.04%. On the other 
hand, the weighted combination of these measure- 
ments gives (|A*n.„|) = 10.0 ± 2.7°, significantly 
different from 0. This finite misalignment probes 
the characteristic timescale of the neutron star 
birth kick, adopting the Spruit & Phinney (1998) 
picture of rotational averaging. As discussed in 
Romani & Ng (2003), the kick timescale con- 
straints are very sensitive to the initial spin peri- 
ods of the individual pulsars. We defer a detailed 
comparison with kick models to a later communi- 
cation. 

A reliable measurement of the spin axis inclina- 
tion ( can also be particularly valuable for these 
young pulsars. Many of these objects are high 
energy (hard X-ray and 7-ray) emitters and the 
modeling and interpretation of the pulse profiles is 
quite sensitive to £ (Romani & Yadigaroglu 1995). 
Radio techniques (polarization sweeps and pulse 
width fitting) are often used to estimate £, but 
these are subject to substantial interpretation un- 
certainties. For example polarization sweep results 
are affected by 90° mode jumps in various pulse 
components, which can often only be resolved with 
single pulse studies. Perhaps this complexity is 
not surprising, since the radio emission, at rela- 
tively low altitude, is sensitive to higher order mul- 
tipoles and the details of the magnetic polar cap 
structure. In the wind zone all such details are 
likely lost. If high energy emission is generated 
in outer magnetosphere (~ 0.1 — 0.3r^c) g a P s , 
then since r^c /r/vs is large, high order multipoles 
should die away and the pulse profiles should be 
sensitive only to magnetic inclination and (. 

We list some radio inclination estimates in Ta- 
ble 2 - however we caution that substantially 
different values are available in the literature 
for many of these pulsars. For Crab and PSR 
J0538+2817 the agreement with our fit £ appears 
good. For Vela and PSR B1706-44 the discrep- 
ancies are large. Interestingly, Helfand, Gotthclf 
& Halpern (2001) match the axis ratio of the 
projected PWN torus by eye and find ( = 55° in 
good agreement with the radio estimate; however 
this value is very strongly excluded in our fits. 



Our relative (s for Vela and PSR B1706-44 do 
however make sense in the outer magnetosphere 
picture (Romani & Yadigaroglu 1995), with PSR 
B 1706-44 at smaller £ producing a narrower dou- 
ble 7-ray pulse and a larger phase delay from the 
radio. We can further make some predictions for 
objects not yet observed at 7-ray energies. If our 
PWN £ estimates hold up, we would expect 7-ray 
emission from PSR J2229+6114 to show a merged 
double pulse, somewhat narrower than that of 
PSR B1706-44 (as appears to be the case in the 
hard X-ray), while PSR J0538+2817 should show 
a wide, Vela-like double pulse. PSR J1930+1852 
in G54. 1+0.3 may be very faint in the 7-ray since 
its small 180° — ( = 33° suggests that any outer- 
magnetosphere 7-ray beams miss the Earth line- 
of-sight. Finally PSR B0656+14 which is many 
times fainter in > lOOMeV 7-rays than expected 
from its spin-down luminosity is widely believed to 
be viewed nearly pole-on. In this orientation the 
strong outer-magnetosphere 7-ray beams would 
not be visible, although we should see 7 emission 
from the pair production fronts in the radio zone 
above the polar cap. New CXO imaging may al- 
low a quantitative fit to the PWN, supporting the 
apparent small £. 

Also physically interesting are the estimates for 
the post-shock velocities (3. For the Crab neb- 
ula, our fit value compares well with the (3 ~ 0.5 
found for the motions of wisps near the torus (Hes- 
ter, et al. 2002), although clearly our very small 
statistical fit errors must under-predict the true 
uncertainty. Moreover, it is puzzling that we get 
slightly larger (3 for the outer ring. One certainly 
expects the flow speeds to drop rapidly in the outer 
nebula and CXO/ HST data do show pattern ve- 
locities as small as (3 ~ 0.03 in the outer parts 
of the X-ray nebula. For Vela, no estimates of (3 
from torus motions have yet been published, but 
Pavlov, et al. (2003) find (3 ~ 0.3-0.6 for features 
in the outer jet, which bracket the torus (3 found 
here. The interpretation of the outer jet speed is 
apparently complicated by a varying orientation 
with respect to the line of sight. Likely the in- 
ner jet /counter- jet provide a cleaner comparison 
with the torus (3; several authors have noted the 
larger counter-jet brightness suggests that it is ap- 
proaching the observer. Using the time-averaged 
image in Figure 1, measuring a 5" length equidis- 
tant from the pulsar on each jet and subtracting 
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the interpolated background from either side, we 
find a counter-jet /jet count ratio Jb = 2.3. For 
a continuous, time-average jet of photon index T 
the the Doppler boosting ratio is 

/ B = [(l+/3cosC)/(l-/3cosC)] r+1 

(the power T + 2 applies for isolated bright blobs). 
Pavlov, et al. (2003) report a inner jet spectral 
index r w 1.1, which with our fit ( gives (3j = 0.45, 
in remarkably good agreement with our torus (3. 
It is worth noting, however, that even these inner 
jets show significant variation, so a time-resolved 
study of pattern speeds and brightness is likely 
needed for a precise jet (3. 

No clear pattern yet emerges from the (3 esti- 
mates presented here, although there is a weak 
anti-correlation between the light cylinder mag- 
netic field and (3 (correlation coefficient « —0.3). 
If significant, this may suggest an anti-correlation 
between pair multiplicity, expected to be large in 
the narrow gaps of high E, short period pulsars, 
and the wind speed. Ultimately, comparison of 
(3 in different PWN components at different lati- 
tudes, and between PWNe, promises to become an 
important new probe of pulsar electrodynamics. 

For Vela, and a few fainter objects, the presence 
of a double ring already suggests some increase in 
pair multiplicity away from latitude 0°. If we de- 
project the ring separation as seen from the pulsar, 
we get a brightness peak (mid-plane of the torus) 
at co-latitude 

<9 Tor = tan- 1 (2rsinC/d) « 74°. 

This is about 10° larger than our best fit (. For 
Vela, the polarization sweep rate maximum sug- 
gests a magnetic axis impact angle (3 = £ — a = 
—4°. Notice that this is on the same side of the 
line-of-sight (a > £) as our fit torus. Thus we 
might plausibly associate the two tori with a near- 
radial outflow of high pair multiplicity plasma 
from the magnetosphcre open zones. Perhaps fur- 
ther observation and modeling could distinguish 
between pairs produced in a vacuum gap near the 
star surface (polar cap) with plasma produced in 
an outer magnetosphcre gap above the null-charge 
surface. For example, the outer magnetosphcre 
gaps should populate field lines at angles > a, con- 
sistent with the observed 9tot > a. If a is close to 
tt/2 (e.g. Crab, PSR J0538+2817) the pair plasma 



from the two poles should merge, leading to a sin- 
gle thicker torus. 

In conclusion, the ubiquitous azimuthally sym- 
metric torus (+jet) structures seen about young 
pulsars provide an important new probe of the 
viewing geometry. We have described a procedure 
for fitting simple geometric models to X-ray im- 
ages that match the gross structure of the central 
regions of many PWNe and provide robust esti- 
mates for model parameters. These models clearly 
do not capture all of the rich, and dynamic, struc- 
tural details seen in the brighter nebulae, such 
as Crab and Vela. However the fitting procedure 
does reduce biases of by-eye 'fits' and does allow 
extraction of geometrical parameters from quite 
faint objects. Accordingly, such fits should allow 
us to measure spin axis orientations for a larger set 
of pulsars and use these angles to probe models of 
the neutron star kick and of pulsar magnetosphcric 
physics. 

This work was supported in part by CXO grant 
G02-3085X. We thank Walter Brisken for discus- 
sion of recent pulsar astrometry. 
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Fig. 2.— Pulsar Wind Tori - C*XO ACIS images 
(left) and best-fit models (right) for Crab, Vela, 
G54.1+0.3, PSR J2229+6114, PSR B1706-44 and 
PSR J0538+2817 (left-right, top-bottom). 
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Table 1 



object 



r(") 







PS/torus sep(") 



0.490±g;ggg -/1.0 x10 s 



Crab (inner) 
Crab (outer) 
Vela 

SNR G54.1+0.3 
PSR J2229+6114 
PSR B1706-44 
PSR J0538+2817 



124.0±8:1 61.3±g;l 15.601°;°^ 3.0* 
126.31±g:g 63.03lg : g 41.331^ 5.9* 
130.63^07 63 - 60 -o:o5 21.25±g;g| 3.0* 



911$ 
103±2 
17512 
155l| 



1471^ 
46l» 
551^ 
9911 



4-6±g:i 

».O_ 2 

o c+0.2 
°-°-0.1 

fi 3+ 10 

°-' J -0.7 



0.74 

i 7+0.3 
i -'-0.7 



U.OOU_ _qoi 
44 +0.004 



i.ilg:i o.62i8:8* 

2.5* 0.49l8;82 



n fi>s+ 003 



•+0.3 n c/i+0.09 
-0.7 u - oq: -0.08 



-/l.l x 10 7 
-/1.3 x 10' 
1701/602 
2221/1113 
384/168 
2442*/52 



f n.6ii8:g 



held fixed in the global fit. 



Table 2 

Proper Motion and Spin Axis Angles 



Pulsar tf |A* n .„| C Cfl 



B0525+21 


124.0 ±0.1 


292 ± 10 


12 ± 10 


61.3 ±0.1 


62° 


B0656+14 




93.1 ±0.4 




small 


16 b 


J0538+2817 


155 ±8 


328± - 4 


7±9 


99 ±8 


97 c 


B0833-45 


130.6 ±0.1 


302 ±4 


8.6 ±4 


63.6 ±0.1 


56 d 


B1706-44 


175 ±4 


160± - 10 


15± - 11 


55 ±0.2 


75 e 


B1951+32 


85± ~ 5 


252 ± 7 


13± - 9 







References: f see text, a see Romani & Yadigaroglu (1995), b Lync & Manchester (1988), c Kramer, et al. 
(2003), d Krishnamohan & Downs (1983), e S. Johnston, priv. comm. 
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